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We study the magneto-optical (MO) response of polar semiconductor BiTel with giant bulk Rashba 
spin splitting at various carrier densities. Despite being non-magnetic, the material is found to yield 
a huge MO activity in the infrared region under moderate magnetic fields (< 3T). By comparison 
with first-principles calculations, we show that such an enhanced MO response is mainly due to the 
intraband transitions between the Rashba-split bulk conduction bands in BiTel, which give rise to 
distinct novel features and systematic doping dependence of the MO spectra. We further predict 
an even more pronounced enhancement in the low-energy MO response and dc Hall effect near the 
crossing (Dirac) point of the conduction bands. 

PACS numbers: 78.20.-c,78.20.Bh,78.20.Ls,71.70.Ej 



The spin-orbit interaction (SOI) plays a crucial role 
in the rapidly evolving field of spintronics ffl-Q. The 
principal importance of SOI is in its ability to intrinsi- 
cally couple the electron spin with its orbital motion, and 
hence produce many novel phenomena such as the spin 
Hall effect [|| and spin Galvanic effect || . In the presence 
of an external magnetic field, SOI can effectively medi- 
ate the interaction between photons and electron spins, 
thereby leading to interesting magneto-optical (MO) ef- 
fects, e.g. non-linear Kerr rotation Wi whereby the polar- 
ization plane of the linearly-polarized light in reflection 
is rotated as a consequence of SOI. 

In practice, most materials under magnetic field ex- 
hibit a rather complicated MO response. This is because 
usually too many energy bands are involved in the op- 
tical excitations, which, given that SOI is also at work, 
leads the respective inter- and intraband optical tran- 
sitions to produce complex spectra. This accordingly 
prevents a comprehensive understanding of the role of 
SOI on the MO response of the given materials using the 
available theoretical models. In contrast, semiconduc- 
tors with Rashba-split conduction/ valence bands appear 
to be ideal systems for studying MO effects, as they have 
a rather simple spin-dependent multi-band scheme. How- 
ever, these systems have rarely been studied up to now, 
mainly because they usually show a very weak Rashba 
spin splitting (RSS) which cannot be resolved experimen- 
tally, and also because RSS had been found only in the 
two-dimensional electron-gas or metallic systems formed 
at the surface or interface where the MO effect is hardly 
detectable. 



This situation has been greatly improved recently by 
the discovery of the giant bulk RSS in the polar semi- 
conductor BiTel. The angle-resolved photocmission spec- 
troscopy (ARPES) || has revealed that the bulk con- 
duction bands in BiTel are subject to a large RSS (see 
Fig. 1(a)), well describable by the 2D Rashba Hamilto- 
nian Hr = p 2 /2m+ Ae^ • (S x p), near the time-reversal 
symmetry point A, where e z is the direction of the po- 
tential gradient breaking the inversion symmetry, and S 
and p are the spin and momentum operators, respec- 
tively. This leads to a Dirac-like band dispersion near 
p = and a splitting increasing linearly with |p|. The 
succeeding first-principles calculations || [l0| and the op- 
tical conductivity spectra JTT| have been consistent with 
the APRES data and further revealed that not only the 
bottom conduction bands (BCBs) but also the top va- 
lence bands (TVBs) are subject to a comparable RSS in 
BiTel. This condition allows several distinct inter- and in- 
traband transitions between these two sets of states. The 
reduced dimensionality together with RSS introduces a 
well-defined singularity in the joint density of states at 
the band edge. Therefore, this material is a promis- 
ing candidate to host the enhanced spin-charge coupling 
and/or the magnctoclcctric effect with the possible ap- 
plications to spintronics. From this viewpoint, the MO 
effect is an intriguing issue in BiTel; the novel interplay 
between RSS and external magnetic field has been ex- 
tensively investigated from theoretical side jl2|-[l5| but 
hardly ever probed experimentally in a real compound. 
Of particular interest are the states near the band cross- 
ing (Dirac) point formed by RSS, as they turn into novel 



2 



Landau levels under magnetic field similar to the case of 
graphenc Thus, it is important to know their im- 

pact on the MO response, e.g. by tuning the Fermi level 
Ep across the Dirac point. In this Letter, we accordingly 
study in detail the MO properties of BiTel both experi- 
mentally and theoretically. It is shown that due to giant 
RSS of bulk bands, the material exhibits a huge enhance- 
ment in its MO response in the infrared spectra under a 
moderately low magnetic field of 3 T. This enhancement 
is further anticipated to be even more pronounced in the 
low-energy region near the Dirac point of BCBs. 

BiTel has a trigonal crystal structure with C% v symme- 
try composed of consecutive Bi, Te, and I layers stacking 
along the c-axis of the crystal || ^| . Due to ionicity and 
covalency of Bi-I and Bi-Te bonds, the bulk material pos- 
sesses an intrinsic polarity along the c-axis which, when 
coupled with the strong atomic SOI of Bi, leads to a 
huge bulk RSS. The calculated band structure of BCBs 
and TVBs around the band gap (along H — A — L di- 
rection), together with the possible inter- and intraband 
transitions, is shown in Fig. 1(a). 
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FIG. 1: (Color online) (a) Calculated electronic band struc- 
ture of BiTel with the a, /3 (intraband), and 7 (interband) 
transitions for a given Fermi level Ef- Ed depicts the position 
of the Dirac point, (b)-(c) carrier density determined from the 
plasma frequency (panel (b), circles) and Hall response (panel 
(c), crosses, see the related discussion in the text). Vertical ar- 
rows indicate the samples investigated experimentally, while 
red, blue, and dashed curves of panel (c) represent the cal- 
culated Hall response with/without taking into account SOI, 
and the conventional relation Rh = — l/(n c e), respectively. 
The inset of panel (b) shows the measured low-temperature 
reflectivity of sample #1. (d) Fermi level as a function of the 
carrier concentration n c . 

The MO effects arc discussed by analyzing the ele- 
ments of the optical conductivity tensor. For the deter- 
mination of the off-diagonal components, first the reflec- 



tivity spectra were measured over a broad energy range 
(E = 10 mcV-40 cV and E = 10 mcV-6 eV at room 
and low temperatures, respectively); for measurements 
above 6 eV, the synchrotron radiation from UVSOR at 
the Institute for Molecular Science, Okazaki, was uti- 
lized. The diagonal part a xx has been derived through 
the Kramers-Kronig transformation. As the second step, 
the complex magneto-optical Kerr angle <S>k was mea- 
sured by a magneto-optical Kerr effect (MOKE) setup: 
The combination of a Fourier-transform infrared spec- 
trometer, a polarization modulation technique jl6j af- 
forded by a ZnSe photoelastic modulator, and a room 
temperature-bore magnet has enabled the direct deriva- 
tion of the off-diagonal conductivity a xy via the following 
relation: 
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The applicable energy range of the setup is 80 — 550 meV. 
The experiments were carried out in reflection configura- 
tion with nearly normal incidence on cleaved surfaces of 
the a6-plane. The low-temperature measurements were 
performed at T = 10 K, while the external magnetic field 
of B = ±3 T has been applied perpendicular to the sam- 
ple surface in the case of the MOKE measurements. 

To analyze the respective MOKE data, the longitu- 
dinal conductivity a xx has been derived from the Kubo 
formula as described in JTT| , while the transverse conduc- 
tivity o~ xy has been calculated to linear order in the ex- 
ternal magnetic field B using the Fukuyama formula [[l7| , 
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where G = [ie n + Ep — H]^ 1 is the one-particle thermal 
Green's function, e„ = e„ + ^sgn(e„) with e„ = (2n + 
l)7rfcsT, r the damping constant, and G_ is obtained 
from G by replacing e„ — >• e n + lui with we = 2^7rfcsT; m 
denotes the electron mass, V the volume of the system, 
and Vi (i = x, y) the velocity operator represented in 
matrix form as 



(n\vi(k)\m) = j^(^ n 



dH(k) 



(3) 



where \n) corresponds to the nth eigenstate of H(k). 
These matrix elements were computed by downfolding 



the ab-initio Hamiltonian |18 to a low-energy tight bind- 



ing model, which was constructed for the 12 valence and 
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6 conduction bands around the band gap M . To evaluate 
the formula (||) , the sum over Matsubara frequencies was 
transformed into a contour integral along the branch cuts 
of the Green's functions |l^] . After performing the ana- 
lytic continuation iuip — > uj, the contour integral as well 
as the k momentum integral were evaluated numerically. 
In our calculations, T is fixed at 12.8 meV (this value has 
already been proven to work well for the analysis of the 
diagonal conductivity spectrum a xx (io) p"l|). 

The carrier densities n c of different samples were first 
estimated from the plasma edge of the free carriers in 
the afr-plane (depicted in the inset of Fig. 1(b)) using 
y / n c e 2 /(m*e n e oc ), where m* is the effective carrier mass, 
and eoo the high-frequency (background) dielectric con- 
stant as determined by the ARPES and previous optical 
measurements, respectively || [y]. The effect of carrier 
doping may be treated within the rigid band approxi- 
mation and then the carrier concentration dependence of 
the Fermi energy shows a monotonic behavior as plotted 
in Fig. 1(d). To study the effect of carrier doping on 
MO response, we targeted the samples with representa- 
tive carrier concentrations, whose Ep are above (sample 
#1 and #2), around (#3), and below (#4) the Dirac 
point. The high-quality single crystal was grown by the 
chemical vapor transport method and n c was controlled 
by the growth conditions for samples #1 and #3. On the 
other hand, the Bridgman method has been employed as 
described in || for samples #2 and #4, and n c has been 
tuned by doping Sb and Cu, respectively. As indicated 
in Fig. 1(b) the corresponding carrier densities estimated 
from the plasma frequencies are n c ~ 7.5, 5.4, 3.4, and 
0.9xl0 19 /cm 3 for samples #1-4, respectively. Another 
way to determine the carrier concentration would be via 
Hall measurements. Importantly, however, in the case 
of BiTel the relation Rh = —l/(n c e) does not hold. 
Figure 1(c) shows the calculated dc Hall response as 
a function of n c with and without taking into account 
SOI, which clearly demonstrates that the conventional 
Hall data analysis under(over)estimates the carrier den- 
sity above (below) the Dirac point. The origin of the 
large enhancement of B / p xy (or the suppression of Rh) 
around the Dirac point is due to the coexistence of elec- 
tron and hole pockets, or equivalently to the donuts-like 
Fermi surfaces in the RSS band. In the case of sam- 
ples #1-3 the results of magnetotransport experiments 
in conjunction with the calculations can give the carrier 
density values n c in accord with the values estimated 
from the plasma edge (Fig. 1(b)). However, this is not 
the case for sample #4, where the complex characteris- 
tics of the Fermi surface makes the simple estimate based 
on the plasma frequency less adequate. Consequently, in 
the following analysis we use n c ~ 0.06 x 10 19 /cm 3 as the 
carrier density of sample #4 estimated from the realistic 
band structure model. 

The experimental spectra for the a6-plane of sample 
#1 (with Ep well above the Dirac point) are plotted in 
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FIG. 2: (Color online) Comparison between the optical con- 
ductivity spectra obtained experimentally (left panels, (a)) 
and theoretically (right panels, (b)) for sample #1. Colored 
lines represent the actual curves, while gray color illustrates 
the Drude-Lorentz response based on the fit to the low energy 
part of the experimental a xx (see the text for details). The 
calculation of a xx (dashed line) does not include the contri- 
bution of free carrier excitations, so this part has been added 
manually by the fit for easy comparison. 



Fig. 2. Our theoretical results, also shown in Fig. 2, 
closely reproduce the experimentally observed spectra 
both in magnitude and energy position, including a xx (ui) 
as well as the real and imaginary parts of a xy (uj). Two 
distinct features in the MO response are discerned in the 
<J xy {uj) spectra around 0.2 eV and 0.4 eV, in addition to 
a Drudc-Lorcntz response dominating the region below 
0.1 eV. These two resonance structures well correspond to 
the intraband transitions a and (3 (see Fig. 1(a)) assigned 
in the a xx (uj) spectra |Q. The observed MO response 
(up to ±5il _1 cm _1 at 0.1-0.5 eV) at 3 T is remarkably 
large for such a non-magnetic system. For a comparison, 
the interband contribution to the transverse conductivity 
for typical nonmagnetic semiconductors InSb, InAs, Ge, 
or GaAs at 3 T is of the order of 0.01-0.3 Sl^cm -1 as cal- 
culated from the Faraday rotation [^(J . Even in the case 
of the ferromagnetic (Gd .95Ca .05)2Mo 2 O7, the MO sig- 
nal in the 0.1-1 eV energy range governed by the Mo Ad 
intraband transitions is an order of magnitude smaller, 
and only the spin chirality induced contribution present 
in Nd2Mo207 is comparable with the present system pl| . 
The observed MO activity at 3 T can well compare with 
the typical magnitude of the MO response {a xy ) coming 
from the p — d charge-transfer excitations in oxide ferro- 
magnets H- In BiTel, because of RSS of BCBs and 
TVBs, the optical transitions between the Rashba-split 
bands are allowed. Since these transitions have different 
dipolc matrix elements for left and right circularly polar- 
ized light under applied magnetic field, a magneto-optical 
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FIG. 3: (Color online) Systematics of spectral features with 
changing the Fermi energy Ef- Ef decreases monotonically 
from sample #1 towards sample #4, crossing the Dirac point 
near sample #3. Vertical broken lines represent the positions 
of the structures identified with the intra- (a and /?) and in- 
terband (7) transitions (see the related discussion). 



response results which scales linearly with the strength 
of SOI [^3j. Hence, the large magneto-optical response 
of BiTcI is a direct consequence of the gigantic bulk RSS 
in this material. This is an astonishing effect of the SOI, 
since the magnetization of BiTcI at B = 3 T is merely 
of the order of 10~ 4 /jb /Bi [Q, several orders of magni- 
tude smaller than the spontaneous magnetization in the 
ferromagnets mentioned above. 

To have a better insight into the consequences of RSS 
on the MO properties, the contribution of free carrier 
excitations has been subtracted by assuming a simple 
Drudc model expressed with use of the cyclotron fre- 
quency uj c = eB/m*, and the relaxation time r — Ti/Y. 
The dc conductivity a xx {u> = 0,B = 0) was determined 
from the transport experiments, then r, the only free 
parameter, was chosen so as to fit the low energy part 
of the experimental diagonal conductivity spectra. The 
Drude-Lorentz curves are also depicted in Fig. 2. Fig- 
ure 3 shows the resulting spectra of a xx and a xy , whose 
Drude-Lorentz components are subtracted with the sim- 
ilar analysis, for all the samples investigated. The values 
of the effective mass (m* = 0.18mn) and the damping 
constant obtained from these fits are in good agreement 
with those found by the ARPES measurement [pi as well 
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FIG. 4: (Color online) Observed energies (Et r ) of intra- and 
interband transitions in BiTel with changing carrier density 
(n c ). Solid lines (symbols) represent the calculated (observed) 
transition energies (energy position of the structures found in 
the MO spectra). Crosses and open circles correspond to the 
features identified with the intra- and interband transitions, 
respectively. 



as those used for the aforementioned theoretical model. 

Next, we explain the evolution of different spectral fea- 
tures found in the MO response with varying carrier den- 
sities. Starting with sample #1, in the a xx (oj) spectra 
one can identify the characteristic features of the intra- 
band transitions as found previously pd| ] , namely the ver- 
tical transitions to (from) the Ep as depicted by a (/3) in 
Fig. 1(a). In external magnetic fields the Zccman split- 
ting is small compared to RSS, nevertheless the matrix 
elements of the intraband transitions become different for 
left and right circularly polarized lights. The dispersive- 
like features found in the real part of a xy (ui) are the con- 
sequence of these processes, giving a clear evidence for 
the strong SOI in BiTcI. Lowering the Fermi level (Ep) 
to the level of sample #2 the a band shifts toward lower 
energy with decreasing the magnitude, while the /3 band 
remains at almost the same position. Further variation 
of Ep (as in samples #3 and #4) causes slight shifts of 
j3 as well, while a is no more disccrnable in the mea- 
sured frequency range. On the other hand, the spectral 
feature characteristic of the interband transitions 7 (see 
Fig. 1(a)) enters into the detectable window from higher 
energies for samples #3 and #4, observed as a sharp rise 
in a xx (uj) at the edge of the interband transition and as a 
dispersive spectral shape in the a xy (oj) spectra. Figure 4 
shows the calculated (solid lines) and observed (symbols) 
results for the relevant transition energies as a function 
of carrier density n c . The agreement between the experi- 
mental results and the theoretical predictions is excellent 
in the case of a and j3 transitions (that are assigned to 
the intraband excitations), and reasonably good in the 
case of the 7 (interband) transitions. 

To summarize, we have studied the magneto-optical 
response of BiTel with a large Rashba spin splitting by 
systematically changing the position of the Fermi level 
around the band-crossing (Dirac) point. Given that 
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BiTcI is a non-magnetic system, the observed MO re- 
sponse arising from the intraband transitions is found to 
be huge compared with conventional (spin degenerate) 
semiconductors. A theoretical model based on the cal- 
culated band structure has been constructed, which can 
quantitatively account for the experimental results, and 
also predict the significant impact of the Dirac point on 
the dc Hall effect. 
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